Evidence for the formation of a Mott state in potassium-intercalated pentacene 
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We investigate electronic transport through pentacene thin-films intercalated with potassium. 
From temperature-dependent conductivity measurements we find that potassium-intercalated pen- 
tacene shows metallic behavior in a broad range of potassium concentrations. Surprisingly, the 
conductivity exhibits a re-entrance into an insulating state when the potassium concentration is 
increased past one atom per molecule. We analyze our observations theoretically by means of 
electronic structure calculations, and we conclude that the phenomenon originates from a Mott 
metal-insulator transition, driven by electron-electron interactions. 



PACS numbers: 73.61.Ph, 71.20.Tx, 71.30.+h 

I. INTRODUCTION 

Pentacene (PEN) is a conjugated molecule very well- 
known in the field of plastic electronics for its use in 
high-mobility organic thin-films transistors^!. Plastic 
electronics applications rely on the fact that at low den- 
sity of charge carriers pentacene films effectively behave 
as weakly doped semiconductors^^. In this regime, 
which is studied extensively, the interactions between 
charge carriers can be neglected. However, the opposite 
regime of high carrier density has remained virtually un- 
explored. Since pentacene froms a molecular solid with 
narrow bandwidth it can be expected that at high den- 
sity the simple assumptions of independent electron-band 
theories break down and the electronic correlations de- 
termine the electronic properties of the materiaP, as it 
happens in other molecular systems. The origin of cor- 
related behavior in these systems is the competition be- 
tween the energy gained by delocalizing the ir electrons 
(given by the electronic bandwidth) and the Coulomb 
repulsion between two carriers on the same molecule. If 
the repulsion energy is larger than the one gained on 
derealization, then the electrons become localized and 
a Mott metal-insulator transition takes place. Among 
the most studied molecular systems in the high carrier 
density regime are the intercalated Ceo crystals-^ and the 
organic charge transfer salts, where the electronic inter- 
actions lead to the appearance of highly correlated mag- 
netic ground states and unconventional superconductiv- 

Our goal is to investigate pentacene compounds at high 
carrier density -of the order of one carrier per molecule- 
and to show that electron correlation effects are crucial 
to understand the resulting electronic properties. To this 
end, we have studied electronic transport through high- 
quality pentacene thin-films similar to those used for the 
fabrication of field-effect transistors. In order to reach 



high densities of charge carriers we intercalate the pen- 
tacene films with potassium atoms, to form K^PEN. Sev- 
eral past experimental studies have addressed the possi- 
bility to chemically dope pentacene thin films through 
the inclusion of alkali atoms and iodine. For all these 
compounds the structural investigations have shown that 
large concentrations of atoms (up to three iodine atoms 
per pentacene molecule) can intercalate in between the 
planes of the pentacene molecular films. Similarly to the 
case of intercalated C^J^l the alkali atoms donate their 
electrons to the Lowest Unoccupied Molecular Orbital 
(LUMO) of the pentacene molecules, whereas the iodine 
donate holes to the Highest Occupied Molecular Orbital 
(HOMO)), enabling th e control of the co nductivity of 
the films. Earlier studie d 11 1 12 ' 13 ' 14 1 15 1 16 1 17 ' 18 1 indicate that 
upon iodine or rubidium intercalation the conductivity 
of pentacene films can become large (in the order of 100 
S/cm) and exhibit a metallic temperature dependence. 
The experiments so far, however, have not led to an un- 
derstanding of the doping dependence of the conductivity 
(e.g. how many electrons can be transferred upon inter- 
calation) or about the microscopic nature of electrical 
conduction of doped pentacene in the high carrier den- 
sity regime. 

In this paper we present the first experimental inves- 
tigation of the evolution of the temperature-dependent 
conductivity of potassium-intercalated pentacene thin 
films, with increasing their doping concentration. We 
find that, upon K intercalation, PEN films become metal- 
lic in a broad range of doping concentrations, up to 
KiPEN, after which the conductivity re-enters an insu- 
lating state. Our experiments also show that the struc- 
tural disorder of PEN films plays an important role on the 
transport properties of K X PEN, as films of poor struc- 
tural quality do not exhibit metallic behavior. The anal- 
ysis of our data shows that at high carrier density the con- 
ductivity of KiPEN cannot be described in terms of indc- 
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pendent electrons filling the molecular band originating 
from the LUMO. Rather, our observations are consistent 
with the formation of a Mott insulating state, driven by 
electron-electron interactions, as we show theoretically 
by calculating the electronic structure of KiPEN. 



ing (RBS). As shown in the top inset of Fig.[T]the ratio of 
K atoms to PEN molecules, Nk/Npen, increases linearly 
with increasing the doping time, as expected. Deviations 
from linearity -approximately 10-20%- are due to inho- 
mogeneity of the potassium concentration. 



II. PREPARATION OF K^PEN FILMS 

Our choice of working with pentacene thin-films (as 
opposed to single crystals) is motivated by both the 
relevance for applications -control of doping in organic 
semiconductors is important for future plastic electronic 
devices- and by the difficulty to grow crystals of alkali- 
intercalated pentacene, which has so far impeded this 
sort of investigations. All technical steps of our exper- 
imental investigations including the deposition of pen- 
tacene films, potassium intercalation, and temperature- 
dependent transport measurements have been carried out 
in ultra-high vacuum (UHV) (10~ n mbar) in a fashion 
that is similar to our previous studies of intercalated ph- 
thalocyanine films^. The use of UHV prevents the oc- 
currence of degradation of the doped films over a period 
of days. 

As with the phthalocyanines, the PEN films (~ 25nm 
thick) were thermally evaporated from a Knudsen cell 
onto a silicon substrate kept at room temperature. In 
order to minimize the parallel conduction through the 
silicon, we use a high resistive silicon-on-insulator (SOI) 
wafer as substrate. The SOI wafer consists of 2/im silicon 
top-layer electrically insulated by 1/im-thick Si02 layer 
from the silicon substrate^. Ti/Au electrodes (lOnm Ti 
and 50nm Au) were deposited ex situ on the SOI sub- 
strates, see Fig. [TJd. After the deposition of the elec- 
trodes and prior to loading the substrate into the UHV 
system, a hydrogen-terminated Si surface was prepared 
by dipping the SOI surface in a hydrofluoric acid solu- 
tion and rinsing in de-ionized water. The use of such a 
H-terminated silicon surface proved necessary to achieve 
sufficient quality in film morphology, as we will discuss 
in the next section in more detail. 

Special care was taken to chemically purify pentacene 
prior to the film deposition. As-purchased pentacene 
powder was purified by means of physical vapor deposi- 
tion in a temperature gradient in the presence of a stream 
of argon gas as described in Ref!^. After this step, the 
pentacene powder was loaded in the Knudsen cell in the 
UHV system and was further purified by heating it at a 
temperature just below the sublimation temperature for 
several days. The film thickness was determined by cal- 
ibrating the pentacene deposition rate ex-situ, using an 
atomic force microscope (AFM). 

Potassium doping was achieved by exposing the films 
to a constant flux of K atoms generated by a current- 
heated getter source. The source was calibrated and the 
potassium concentration determined by means of an el- 
emental analysis performed on PEN films doped at sev- 
eral doping levels, using ex-situ Rutherford backscatter- 



III. TRANSPORT PROPERTIES OF K X PEN 

A. Electronic transport through high 
structural-quality K^PEN films 

The conductance of K^PEN films is measured in situ in 
a two terminal measurement configuration with a contact 
separation of approximately 175 /im (see Fig. [lja). The 
dependence of the conductivity on the potassium con- 
centration, hereafter referred to as the "doping curve", is 
determined for different PEN films, as a function of the 
ratio of K atoms to PEN molecules. The doping curves 
for different samples are very similar, as shown in Fig. [TJi. 
Upon doping, the conductivity initially increases rapidly 
up to a value of a ~ 100 S/cm -in the same range as 
the conductivity of metallic K 3 C 6[ pll. Upon doping fur- 
ther, the conductivity continues to increase more slowly, 
reaches a maximum at a concentration of one K per PEN 
and then drops sharply back to the value of the undoped 
PEN film. All of the more than 40 films that we have 
investigated exhibit a similar behavior. 

The observed suppression of the conductivity of pen- 
tacene films at high doping (for potassium concentrations 
higher than 1 K/PEN) allows us to exclude the possibil- 
ity that the conduction of the intercalated films observed 
in the experiments is due to an experimental artifact, 
for instance the formation of a potassium layer on top 
of the pentacene film. In fact, at doping higher than 1 
K/PEN the measured conductance, and its temperature 
dependence, is essentially identical to what is measured 
for pristine films. 

To understand the nature of conduction of pentacene 
films at high carrier density we measured the tempera- 
ture dependence of the conductivity for different values 
of potassium concentration (see Fig. [2^,) . Pristine PEN 
films have a very low conductivity and the measured 
conductance of undoped films is dominated by trans- 
port through the substrate's 2 /im thick Si top-layer. 
The measured conductivity decreases rapidly with lower- 
ing temperature, as expected, confirming that undoped 
(x—0) pentacene films are insulating. On the contrary, in 
the highly conductive state -for x between 0.1 and 1 -the 
conductance of the films remains high down to the lowest 
temperature reached in the experiments (~5K), indicat- 
ing a metallic state. When the potassium concentration 
is increased beyond approximately 1 K/PEN, the conduc- 
tivity again decreases rapidly with lowering temperature, 
indicating a re-entrance into an insulating state. The 
metallic and insulating nature of pentacene thin-films 
at different potassium concentrations is confirmed by 
measurements of volt-amperometric characteristics [I-V 
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FIG. 1: (a) Conductivity a and square conductance Gu, of 
three different K-doped PEN films, as function of the ratio 
Nk/Npen; under the curves a pentacene molecule. Inset: 
Nk/Npen as a function of doping time. Schematic view (b) 
of our set up and (c) atomic force microscopy image of a high 
quality undoped PEN film showing large crystalline grains. 
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FIG. 2: Temperature dependence of the conductance of potas- 
sium intercalated pentacene films, (a) The colored dots in the 
inset of (a) indicate the doping level at which the tempera- 
ture dependent conductivity measurements with correspond- 
ing color were performed. In black the temperature dependent 
conductivity of the Si substrate is shown. The low tempera- 
ture (5K) I-V characteristics of K^PEN in the (b) conducting 
and (c) highly-doped insulating state. 



curves) at 5 K. For x between 0.1 and 1 the films ex- 
hibit linear I-V characteristics (Fig. [2b, as expected for 
a metal. In the highly doped regime (for x > 1), on the 
contrary, the insulating state manifests itself in strongly 
non-linear I-V curves and virtually no current flowing at 
low bias (Fig. |2p). Therefore, the data clearly show that 
pentacene films undergo a metal-insulator transition as 
the density of potassium is increased past one atom per 
molecule. Since in the overdoped regime the conduction 
occurs through the Si layer of the SOI substrate, it is not 
possible to gain specific information about the properties 
of the insulating KPEN films -for instance, to determine 
the electronic gap from measurements of the activation 
energy of the conductivity- by studying dc transport on 
our samples. 



B. The effect of structural disorder on the 
transport properties of K T PEN 

The high structural quality of the films proves to be 
the essential ingredient necessary to obtain K^PEN films 
which exhibited metallic conductivity. We find that the 
quality of pentacene thin films is highly sensitive to the 
choice of the substrate material and sufficient quality can 



be achieved by using a hydrogen-terminated Si surface. 
To illustrate this important technical point, we show here 
that the structural quality of films deposited on a Si02 
surface has a very large impact on their electronic trans- 
port, with low quality resulting in considerably poorer 
electrical properties. 

Fig. [3^ shows the doping curve of PEN films deposited 
onto 300nm SiC>2 that was thermally grown on a Si sub- 
strate. For these films, the maximum conductivity that 
we measured experimentally is several orders of magni- 
tude lower than the conductivity measured for films de- 
posited on a Si surface. In addition, (on Si02) the con- 
ductivity was always observed to decrease rapidly with 
lowering temperature, i.e. the potassium intercalated 
films are always insulating (see Fig.j3j)). Both the magni- 
tude and the temperature dependence of the conductivity 
that we measured on SiC>2 substrates are comparable to 
results obtained in earlier work reported in the literature. 

We attribute the difference in the electrical behavior 
observed for films deposited on Si and Si0 2 substrates 
to the difference in film morphology, which we have ana- 
lyzed using an atomic force microscope. Figure [4] shows 
AFM images of two pentacene films of similar thickness 
deposited on the hydrogen terminated Si (Fig. and b) 
and on the SiC>2 surface (Fig. [it and d). It is appar- 
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FIG. 3: (a) Conductivity, a measured at room temperature 
as a function of doping time for a 25nm thick pentacene film 
deposited on Si02- (b) Temperature dependence of the con- 
ductivity for a pentacene film grown on Si02 and doped into 
the highest conductivity state. The conductivity is rapidly de- 
creasing with lowering the temperature as it is typical for an 
insulator, (c and d) AFM images of pentacene films grown on 
Si02 and on H-terminated Si. Small and randomly oriented 
grains with large height fluctuations are observed when the 
PEN films are deposited on Si02 (c), whereas PEN films of 
similar thickness deposited on Si consist of large crystalline 
grains with a common relative orientation and only relatively 
small fluctuations in height (d). 
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FIG. 4: (a) Temperature dependence of the resistivity of a 
high-quality K^Ceo film grown on Si. As expected, the re- 
sistivity exhibits a superconducting transition at 18 K. The 
inset shows the doping dependence of the conductivity. The 
conductance peak is typical of K^Ceo- (b) Temperature de- 
pendence of the resistivity of a K^Ceo film grown on Si02 
showing insulating behavior. 



ent that very different morphologies are observed for the 
two substrates. PEN films deposited on Si surfaces ex- 
hibit large crystalline grains with a common relative ori- 
entation and only relatively small fluctuations in height. 
On Si02, on the contrary, the grains are much smaller, 
randomly oriented and they exhibit much larger height 
fluctuations. 

This conclusion is consistent with past studies^ show- 



ing that the growth and morphology of pentacene films is 
strongly influenced by the substrate surface. Specifically, 
for pentacene films grown on SiC>2 , a high density of nu- 
cleation centers was observed, leading to the growth of 
small islands and to a high concentration of grain bound- 
aries. On the hydrogen terminated silicon surface, on 
the other hand, the much smaller density of nucleation 
centers results in significantly larger islands and in a re- 
duced density of grain boundaries. Note that the crit- 
ical influence of the film morphology on the electrical 
characteristics of electron-doped pentacene films is also 
supported by recent experiments studying the conduc- 
tion of rubidium-intercalated pentacene films deposited 
on glass^. In that work, as-doped films exhibited an 
insulating temperature dependence of the conductivity. 
However, by performing a high-temperature annealing 
on the doped films, which results in an improved mor- 
phological quality, metallic behavior was also observed. 

The sensitivity of the morphology of pentacene films to 
the substrate, together with the resulting effects on the 
electronic properties, is common to films of many conju- 
gated molecules. In fact, a similar sensitivity was found 
in our earlier work on the electronic properties of alkali 
doped Metal-Phthalocyanine (MPc) filmgSS. Specifically, 
for films of CuPc, NiPc, ZnPc, FcPc, and MnPc, the 
maximum conductivity which can be achieved upon alkali 
doping when the films are deposited on Si02 substrates 
is several orders of magnitude lower than the conductiv- 
ity measured for films deposited on a Si surface, and has 
always an insulating temperature dependence. Also for 
alkali doped Cgo films, it has been observed that the sur- 
face termination of the substrate affects the morphology 
and the electronic transport properties of the films. As 
illustrated in Fig. [4^, the resistivity of K3C60 films grown 
on Si shows a low resistivity at low temperature and a 
transition to a superconducting state. On the contrary, 
the K3C60 films grown on Si02 have significantly higher 
resistivity, exhibiting thermally activated transport (see 
Fig. HJ)), without a superconducting transition. 



IV. INTERPRETATION IN TERMS OF A 
MOTT STATE OF KiPEN 

The most striking aspect of our observations, namely 
a sharp decrease in the conductivity starting at a car- 
rier concentration of one electron per molecule concomi- 
tant with the re-entrance into an insulating state, has 
not been reported in earlier experiments on intercalated 
pentacene (in which the density of intercalants could not 
be determined 14 ! 15 * 16 ! 17 * 18 !), nor in studies of pentacene 
field-effect transistors with gate electrolytes (in which a 
metallic state has not been observed^. It implies that, 
contrary to the case of pentacene devices used in plastic 
electronics, the electronic properties of pentacene films at 
high carrier density cannot be described in terms of non- 
interacting electrons. In fact, even though it is known 
that pentacene molecules can accept only one electron 



5 



and that d oubly negatively charged pentacene ions do not 
exist! 23 * 24 -! (i.e., in our films charge transfer from the potas- 
sium atoms saturates at 1 K/PEN), a carrier concentra- 
tion of one electron per molecule corresponds to a half- 
filled band, and, for non-interacting electrons, should re- 
sult in a metallic state. Therefore, interactions need to 
be invoked in order to explain our observations. 

An established scenario for the formation of an insu- 
lating state at half-filling is the one of a Mott insu lator 
emerging from strong electron-electron interactions^ 2 ^. 
In a Mott insulator a strong Coulomb repulsion prevents 
two electrons to occupy the same pentacene molecule. 
Since at half-filling the motion of electrons necessarily 
requires double occupation of molecular sites, electron 
transport is suppressed and the system becomes insulat- 
ing. This scenario is usually modelled theoretically using 
a Mott-Hubbard Hamiltonian, which includes a kinetic 
energy term described within a tight-binding scheme, and 
an on-site repulsion term. The Mott state occurs when 
this repulsion (U) is larger than the bandwidth W (de- 
termined by the tight-binding hopping amplitudes t). In 
this case the half-filled band splits into a lower (com- 
pletely filled) and an upper (completely empty) Hubbard 
band, separated by a (Mott) gap of the order of U, when 
the interactions are strong. It is realistic that this sce- 
nario is realized in a molecular solid such as pentacene, 
in which the bandwidth is expected to be small owing to 
the absence of covalent bonds between the molecules. 

To substantiate the Mott-insulator hypothesis we have 
analyzed the electronic structure using density functional 
theory (DFT) calculations to extract the parameters of 
the Mott-Hubbard model. A main difficulty in doing 
this is that the structural knowledge of the intercalated 
films is incomplete, as our ultra-high vacuum set-up is not 
equipped to perform in situ structural characterization, 
and ex situ characterization is impeded by oxidation of 
potassium when the sample is extracted from the vacuum 
system where the films are prepared. Therefore, for the 
DFT calculations we take advantage of the existing struc- 
tural information on intercalated pentacene compounds 
and we determine the stable crystal structure of KiPEN 
using a computational relaxation procedure which refines 
the positions of all the atoms in the unit cell. 



A. Structural Details of KiPEN 

It is well known from previous structural studies 
on pentacene films that the herringbone arrangement 
of the molecules is p reserved when pentacene is in- 
tercal ated with i odine ^ 12 * 13 ! or with different alkali 
atoms^ 15 * 16 ! 17 * 18 !, and that intercalation takes place be- 
tween the pentacene layers. It is also known that inter- 
calation is accompanied by a considerable expansion of 
the unit cell c-axis (by an amount close to the radii of the 
intercalated ions) while the in-plane lattice parameters a 
and b are only minorly affected 1 ^ 1 ^!. 

A reliable estimate of the length of the expanded c- 



axis is given by the sum of the radius of the alkali ion 
and the pristine c-axis parameter: the c-axis lattice con- 
stants that are obtained in this way for, e.g., RbPEN 
and CsPEN are within 2% of the experimental values. 
For KiPEN we construct the lattice parameters starting 
from two different polymorphs (one with c = 14.33 Aand 
the other with 14.53 A), using a K + ionic radius of 1.33 
A. The structures of the two polymorphs are taken from 
the experimental results of Ref.™ They differ slightly in 
the packing of the pentacene molecules, which enables us 
to study the influence of realistic variations of the packing 
on the electronic structure. The actual size of the c-axis 
in the KiPEN is not critical. Our relaxation and band- 
structure computations were checked for values up to 8% 
larger and smaller than the estimated c-axis parameters. 
We found that even such relatively large variations of 
c do not affect our main results (i.e., the values of the 
calculated bandwidth W and on site repulsion U), be- 
cause the interaction between adjacent pentacene layers 
is weak. 

After constructing the unit cell of potassium- 
intercalated pentacene, using the information above to fix 
a, b and c, we refine the positions of the atoms by a com- 
putational relaxation procedure. For all the electronic 
structure calculations we used the Vienna Ab Initio sim- 
ulation package (VASP^P^U w ith projector au gm ented 
waves (PAW)Pand the PW91 density functional. The 
self-consistent calculations were carried out with an in- 
tegration of the Brillouin using the Monckhorst-Pack 
scheme with a 6x6x4 fc-points grid and a smearing pa- 
rameter of 0.01 eV, and a plane-waves basis set with a 
cutoff energy of 550 eV. To determine the stable struc- 
ture of KiPEN all the atoms positions in the unit cell 
are relaxed using a conjugate-gradient method. To avoid 
possible energy barriers we used a number of different 
initial configurations. In the relaxation procedure first 
the forces on the K-ions are calculated and then the K 
positions are relaxed. We observe that the dopants move 
into high symmetry positions in the plane between the 
pentacene layers. In the next step the positions of all 
atoms in the unit cell are relaxed -including the ones of 
the two PEN molecules. The final stable structure is the 
same for all different initial configurations.^ 

The optimized structure of KiPEN is shown in Fig. [5] 
We checked the reliability of the relaxation procedure 
on undoped pentacene and found that the calculated 
structure indeed corresponds to the actual, known crys- 
tal structure of the material. In KiPEN there are two 
inequivalent PEN molecules per unit cell, just as in the 
undoped compound. Intercalation changes the detailed 
molecular orientations in the unit cell, see Fig. [5j but we 
do not observe the formation of superstructures such as 
for instance molecular dimers. For the two distinct pen- 
tacene polymorphs^ for which we have performed the 
relaxation procedure, we found that the conclusions on 
electronic bandwidths and Coulomb interactions that will 
be presented hereafter hold equally well. 
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FIG. 5: Crystal structure of potassium-intercalated pentacene 
KiPEN obtained from ab initio computational relaxation 
with (a) showing the herringbone of the PEN molecules and 
K atoms in the unit cell and (b) is a side view of the stacked 
layers of PEN and K, illustrating the potassium intercalation 
in between the molecular planes. The unit cell parameters are 
a,b,c = 6.239, 7.636, 15.682 A, a,p,j = 76.98, 88.14, 84.42° 



B. Electronic Structure and Electronic 
Correlations in KiPEN 

Figure [6^, shows the DFT band structure of KiPEN 
together with the projected density of states on the pen- 
tacene and potassium orbitals for the polymorph associ- 
ated with the relaxed structure of Fig. [5] The Fermi en- 
ergy lies in the middle of a half-filled band that is entirely 
of pentacene character, originating from its LUMO. The 
potassium derived electronic states are present only at 
much higher energy, demonstrating that little hybridiza- 
tion takes place and that the role of the potassium atoms 
is limited to transferring its electrons to the pentacene 
molecules. The total bandwidth is W=0.7 eV. From a 
tight-binding fit of the band dispersion (see Fig. |6|d) we 
extract the hopping amplitudes tij that enter the kinetic 
energy part of the Mott-Hubbard Hamiltonian 
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FIG. 6: Results of electronic structure calculations for 
KiPEN. (a) Single-particle band-structure 1 ^, with the Fermi- 
level Ef (green line) as the zero of energy. Valence and 
conduction band are indicated by the thick blue lines, (b) 
Carbon (blue) and potassium (orange) projected density of 
states, (c) Tight-binding fit to the valence and conduction 
band (blue thin lines) and the resulting lower and upper Hub- 
bard band from a mean field analysis of the corresponding 
Hubbard Hamiltonian with U = 1.45 eV (red thick lines). 
The arrows indicate the opening of the Hubbard gap. (d) Re- 
ciprocal lattice vectors and the first Brillouin zone of KiPEN. 
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TABLE I: Tight-binding fit parameters to the ab initio 
band-structure of the half-filled conduction/ valence band of 
KPEN . The on-site energy difference is denoted by e and 
the hopping parameters along the a, b and c axes are denoted 
by t. 



(ij) 



where we have two molecules in the unit cell, with on-site 
energy e^, the electron creation (annihilation) operators 
on site i are c| CT {c i CT ), with a the electron spin, h.c. the 

hermitian conjugate, i%i jC — cj a Ci^i n i — Yla n h^ an d 
U the effective Coulomb interaction between two elec- 
trons on the same molecule. The hopping integrals ij* 
are different in different directions and between nearest 
and next nearest neighbor molecules, see Table [i] The re- 
sulting electronic band-structure (Fig. [6^,) displays only 
very minor differences for the two stable polymorphs and 
within our accuracy the tight-binding parameters are the 
same. 

In order to determine the relative strength of electronic 
correlations and to compute the magnetic exchange in- 
teractions, the on-site Coulomb interaction Ubare for two 



electrons on the same pentacene molecule is determined, 
using the techniques described in RefP. For this the to- 
tal energy of neutral and charged pentacene molecules 
is calculated by density functional calculations in the lo- 
cal density approximation using GAMESS with a DZVP 
basis-set 37 . The bare value of the Coulomb interaction 
is found to be Ubare = 3.50 eV. In t he solid this value 
is screened, leading to a lower value £y l 38 l 39 l 40 l. From the 
eigenvalues of the charged molecule that is placed in- 
side a cavity of an homogeneous dielectric medium with 
dielectric constant of 3.3 using the SS(V)PE modePSl 
one finds U — 1.45 eV. We have also performed an inde- 
pendent estimate for the value of U by considering the 
difference between the band gap of pristine pentacene 
from density functional calculations (0.7 eV) and its ex- 
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perimental value (2.2 eV^J), which gives U ~ 1.5. These 
two values, determined in two very different ways, are 
remarkably close. Very similar values for U are found 
also for the second polymorph used in our calculations, 
indicating that these values are not very sensitive to dif- 
ferences in the structure 3 ^!. 

From a straightforward self-consistent mean-field de- 
coupling computation on the resulting Hubbard Hamil- 
tonian the ground state is found to be a Neel ordered 
antiferromagnet, with a charge gap of 1.23 eV. The an- 
tiferromagnetic exchange between neighboring molecules 
in the plane is J = 4t? a±6 w 2 /?7 ~ 290K. This value is 
actually an underestimation of the Heisenberg exchange, 
as certainly a nearest neighbor Coulomb interaction V is 
also present, which has the effect of increasing the value 
of exchange by a factor to U/ (U — V)W We find that 
the coupling between molecules in neighboring planes, 
J± , is four orders of magnitude lower than the in-plane 
J. Consequently KiPEN is a quasi two-dimensional an- 
tiferromagnet. Finally, an antiferromagnetic exchange 
of ~ 40K is also present between in-plane next nearest 
neighbor molecules along the a-axis, leading to a weak 
frustration of the magnetic Neel ordering. 

An alternative procedure to take the electronic corre- 
lations into account, used in band structure calculations 
of strongly correlated transition metal oxides, is to in- 
corporate the atomic Hubbard U directly in a LDA+U 
scheme. For molecular crystals, however, this approach 
would require a novel computational scheme that uses 
a basis set of orthonormal molecular orbitals so that a 
molecular U can be incorporated. 

V. DISCUSSION AND CONCLUSIONS 

Using the results of the electronic structure calcula- 
tions we are now in a position to validate the Mott-state 
hypothesis. With a ratio U/W ~ 2.1, electron-electron 
interactions cause the splitting of the LUMO band and 
the opening of a Mott gap, as shown in Fig. [6j:. The 
gap explains the observed re-entrance to the insulating 
state. Note that the Mott-state scenario also explains 
why the insulating state is only observed for a potas- 
sium concentration of 1.1-1.2 atoms per molecule (and 
not at exactly one; see Fig. [TJi). In fact, at exactly one 
potassium per pentacene molecule, non-uniformity in the 
potassium concentration -estimated to be approximately 
10-20% in our films - effectively dopes the Mott insulator 
causing the conductivity to remain large. However, even 
in the presence of non-uniformity, a potassium concentra- 
tion slightly larger than 1 K/PEN results in a uniform 
electron concentration exactly equal to one electron per 



molecule, since, as we mentioned earlier, only one elec- 
tron (a nd no t two) can be donated to each pentacene 
molecule 1 23 * 24 * It should be noted that imperfections in 
the material, cither due to disorder of the dopants or in 
the molecular arrangements will lead to the presence of 
both disorder in the bandwidth and the presence of a 
disorder potential. In general the physics of disordered 
Mott-Hubbard systems is very ricrF 3 -', but it is not a priori 
clear how relevant disorder will be in the present situa- 
tion as we find that the Mott state in KiPEN is stabilized 
by a substantial electronic gap. 

Since the coupling between pentacene molecules in dif- 
ferent layers is very small and the electron-electron inter- 
action is sufficiently large, the low-energy effective elec- 
tronic Hamiltonian of the K^PEN reduces to the well- 
known two-dimensional iJ-modeP' with t/J w 3 — 4. 
Interestingly, the same t J-model in the same coupling 
regime describes another important class of materials, 
namely strongly correlated cuprate superconductors such 
as L^-zSr^CuCj. An apparent difference between these 
classes of materials is that it is generally excepted that 
in doped organics the formation of lattice-polarons plays 
a very important role. 

We conclude that temperature-dependent transport 
measurements and theoretical calculations consistently 
indicate that at a doping concentration of one potassium 
ion per molecule potassium-intercalated pentacene is 
a strongly correlated Mott insulator, whose electronic 
properties are dominated by electron-electron interac- 
tions. An immediate consequence is the emergence of 
magnetism. Our calculations show that the magnetic 
interactions are dominated by a large, positive magnetic 
exchange J = 4i 2 /[/ ~ 290 K between electrons on 
nearest neighbor molecules in the same pentacene layer. 
We predict that KiPEN is therefore an antiferromagnet. 
In fact, experimental indications for the presence of 
antifcrromagnetism in intercalated pentacene have 
been reported in magnetic susceptibility measurements 
performed in the pastp2, albeit at very low temperature. 



VI. ACKNOWLEDGEMENTS 

This work was supported by the Foundation for 
Fundamental Research on Matter (FOM), the Royal 
Dutch Academy of Sciences, the NWO Vernieuwingsim- 
puls, NanoNed and the Stichting Nationalc Computer- 
faciliteitcn. We are grateful to the FOM Institute for 
Atomic and Molecular Physics (AMOLF) for the RBS 
analysis of our samples. 



1 S. F. Nelson et al, Appl. Phys. Lett. 72, 1854 (1998). Mater. 14, 99 (2002). 

2 C. D. Dimitrakopoulos and P. R. L. Malenfant, Adv. 3 H. Sirringhaus et al, Science 290, 2123 (2000). 



8 



4 G. H. Gelinck et al, Nat. Mater. 3, 106 (2004). 

5 H. Klauk et al, Nature 445, 745 (2007). 

6 G. Brocks, J. van den Brink, A. F. Morpurgo, Phys. Rev. 
Lett. 93, 1464051 (2004). 

7 T. Takenobu et al, Phys. Rev. Lett. 85, 381 (2000). 

8 T. Ishiguro, K. Yamaji, G. Saito, Organic Superconductors 
(Springer- Verlag, Berlin) (1998). 

9 M. Imada, A. Fujimori, Y. Tokura, Rev. Mod. Phys. TO, 
1039 (1998). 

10 L. Forro and L. Mihaly, Rep. Prog. Phys. 64, 649 (2001). 

11 T. Minakata, I. Nagoya, M. Ozaki, J. Appl. Phys. 69, 7354 

(1991) . 

12 T. Minakata, H. Imai, M. Ozaki, J. Appl. Phys. 72, 4178 

(1992) . 

13 T. Ito et al, J. Phys. Chem. Sol. 65, 609 (2004). 

14 T. Minakata, M. Ozaki, H. Imai, J. Appl. Phys. 74, 1079 

(1993) . 

15 Y. Matsuo, S. Sasaki, S. Ikehata, Phys. Lett. A 321, 62 

(2004) . 

16 Y. Matsuo et al, J. Phys. Chem. Sol. 65, 619 (2004). 

17 Y. Kaneko et al, Synth. Met. 154, 177 (2005). 

18 B. Fang, H. Zhou, I. Honma, Appl. Phys. Lett. 86, 2619091 

(2005) . 

19 M. F. Craciun et al, Adv. Materials 18, 320 (2006). 

20 R. W. I. de Boer, M. E. Gershenson, A. F. Morpurgo and 
V. Podzorov, Phys. Stat. Sol. (a) 201, 1302 (2004). 

21 T. T. M. Palstra et al, Phys. Rev. Lett. 68, 1054 (1992). 

22 M. J. Panzer and C. D. Frisbie, J. Am. Chem. Soc. 127, 
6960 (2005). 

23 J. Szczepanski, C. Wehlburg, M. Vala, Chem. Phys. Lett. 
232, 221 (1995). 

24 T. M. Halasinski et al, Phys. Chem. A 104, 7484 (2000). 

25 R. W. Lof et al., Phys. Rev. Lett. 68, 3924 (1992). 

26 G. Kresse and J. Hafner, Phys. Rev. B 47, (R)558 (1993). 

27 G. Kresse and J. Furthmiiller, Phys. Rev. B 54, 11169 
(1996). 

28 G. Kresse and D. Joubert, Phys. Rev. B 59, 1758 (1999). 

29 J. P. Perdew et al, Phys. Rev. B 46, 6671 (1992). 

30 C. C. Mattheus et al, Acta Cryst. C 57, 939 (2001). 
The points in the Brillouin zone are F = (0, 0, 0),X 



(!,o,o),m = (|,|,o),y 

(0, !,§),#= i) andO = (i 0,|) 



(0,|, 0),Z = (0,0,|), N 



' 2 ' 2 / ' V 2 ' 2 ' 2 / w V 2 ' ' 2 f 

32 In an alternative crystal structure for intercalated pen- 
tacene, see RefP^, a still smaller value of the bandwidth is 
found, implying an even larger correlation gap. 

33 A. Hansson, J. Bohlin, S. Stafstrom, Phys. Rev. B 73, 
184114 (2006). 

34 T. Mori and S. Ikehata, J. Appl. Phys. 82, 5670 (1997). 

35 R. Eder, J. van den Brink, G. A. Sawatzky, Phys. Rev. B 
54, R732 (1996). 

36 R. Ruiz, et al, Phys. Rev. B 67, 1254061 (2003). 

37 M. W. Schmidt et al, J. Comput. Chem 14, 1347-1363 
(1993). 

38 G. Giovannetii, G. Brocks and J. van den Brink, Phys. 
Rev. B 77, 035133 (2008). 

39 J. van den Brink, et al, Phys. Rev. Lett. 75, 4658 (1995), 
R. Eder, J. van den Brink and G.A. Sawatzky, Phys. Rev. 
B. 54, R732 (1995). 

40 M.B.J. Meinders et al. Phys. Rev. B. 52, 2484 (1995), J. 
van den Brink, R. Eder and G.A. Sawatzky, Europhys. 
Lett. 37, 471 (1997). 

41 D. M. Chipman, J. Chem. Phys. 112, 5558 (2000). 

42 E. A. Silinsh, V. A. Kolesnikov, I. J. Muzikante and D. R. 



Balode, Phys. Status Solidi B 113, 379 (1982). 

43 D. Belitz and T. R. Kirkpatrick, Rev. Mod. Phys. 66, 261 
(1994). 

44 The data file with the resulting structure for both poly- 
morphs is available directly from the authors. 



